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Abstract. We extend the exact periodic Bethe Ansatz solution for one- 
dimensional bosons and fermions with 5-interaction and arbitrary internal degrees 
of freedom to the case of hard wall boundary conditions. Wc give an analysis of 
the ground state properties of fcrmionic systems with two internal degrees of 
freedom, including expansions of the ground state energy in the weak and strong 
coupling limits. 
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1. Introduction 

One-dimensional quantum gases with two-particle (5-interaction have long been of 
fascination. The most simple model of (^-interacting spinless bosons in a periodic 
box was solved in terms of the Bethe Ansatz by Lieb and Liniger [Q. This quantum 
mechanical model is not only one of the oldest integrable models after the Heisenberg 
spin chain, but arguably also one of the most important test beds for exploring new 
ideas and methods, e.g., the Thermodynamic Bethe Ansatz [2] was pioneered for this 
model. Earlier Girardeau [3] discussed a mapping from strongly repulsive interacting 
bosons to fermions, corresponding to the limit c^oo, where c features as the arbitrary 
interaction strength in the Lieb-Linigcr model. Later in seminal work McGuire 0] 
discussed (5-interaction particles via an optical analogue. Gaudin [Fj and Yang [S] 
then considered spin-i fermions with periodic boundary conditions, the first model 
with internal states. Sutherland [7] applied the nested Bethe Ansatz, which allowed 
the treatment of periodic quantum gases with arbitrary spin by repeated application 
of the Bethe Ansatz, reducing the number of internal states in each step. Gaudin [5] 
solved the model of spinless bosons with hard wall boundary conditions. 

The special form of (5-intcraction is at the heart of the intcgrability of the quantum 
gases. Models tweaking the type of interaction have been considered, but so far are less 
prominent UJ. Non- integrable models, like the harmonically trapped gas with tunable 
interaction strength, have to be treated by approximate methods and simulations. 
One integrable model, the mixture of fermionic and bosonic particles, solved by Lai 
and Yang |10j has long been dormant in the literature, but is enjoying new interest 
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11 in light of the novel recent experiments on low-dimensional quantum gases. Apart 
from isolated work |12| this is the only integrable approach to mixtures of bosons and 
fermions. 

Generalizations of the original models and the evaluation of their physical 
properties have been discussed extensively in several books El Q5I E3 QZI El 
and review articles |19| , mostly concentrating on periodic spinless bosons and periodic 
spin-i fermions. Of direct relevance here are a series of papers examining periodic 
su(2) |201 and su(3) [21] bosons as well as hard wall su(2) bosons [221 via the Bcthe 
Ansatz. Hard wall fermions with two internal states have been treated numerically 
for up to three particles [221, but no Bethe Ansatz solution was given. A technically 
involved model in which the hard walls are of finite height, resulting in tunnelling 
outside the box, has also been considered [23- However, only the case of one particle 
outside the box could be treated and no physical properties where evaluated due to 
the complicated nature of the solution. In 1985 Woynarovich [2H] discussed a hard 
wall spin-i fermion model with an additional internal potential. The Bethe Ansatz 
solution he obtained differs from the standard model by one additional term The 
su(N) pure fermion model with <5-potentials of variable strength at the boundaries, 
acting like transparent walls, was also solved [2EJ by extending this idea. This system 
reduces to the hard wall model considered here for |7o|, \ \ ~ * °°, where 70 and 7^ act 
as boundary potentials. Finally we note that the Hubbard model with open ends was 
also solved exactly \27\ . In an appropriate continuum limit |28| the resulting integral 
equations reduce to the integral equations for the ground state energy of a hard wall 
fermion gas with spin-^. 

The content and arrangement of this paper is as follows. In Section [21 we discuss 
the most general one-dimensional integrable (^-interaction model and its exact solution, 
both for periodic and hard wall boundary conditions. The treatment for periodic 
boundary conditions serves to highlight the key differences in the solution for hard 
wall boundary conditions. In Sections[31and01we give an approximate treatment of the 
Bcthe Ansatz solution for the special case of the spin-i fermionic gas with periodic and 
hard wall boundary conditions via methods developed in References (2H1 HUH EH] ■ The 
ground state energy in the thermodynamic limit is considered via integral equations 
in Section [5] Concluding remarks are given in Section [HI The detailed derivation of 
the coordinate Bethe Ansatz for the first level and the Algebraic Bethe Ansatz for the 
nested structure for bosons and fermions with arbitrary spin and either periodic or 
hard wall boundary conditions is given in the Appendices. 

2. Exact solution 

As already remarked, special cases of the exact solutions presented here have been 
examined for more than 40 years. The one-dimensional ^-interaction quantum gas 
with N particles is described by wave functions . . . , xjv, 01, . . . , on) which are 

eigenfunctions of the Hamiltonian H = 7i © ld\r 1 © ■ ■ ■ © Idvw > where 



Here the space coordinates, Xi € [0, L], i = 1, . . . , N, occur only in the space part TL 1 
which operates on the standard space of square integrable functions L? [0, L] , as usually 
used in quantum mechanics. The spin coordinates Oi take values in Oi € {1, . . . , N'} 
and <Tj acts trivially on the AT'-dimensional vector space Vj. 
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Here we initially solve the problem for an abstract model with N' internal states. 
The name su(N') model results from the symmetry of the underlying integrable model 
of the su(N') chain, but one can make the case that this quantum gas Hamiltonian 
also possesses su(N') symmetry itself. For physical application these internal states 
are usually interpreted as various distinguishing physical properties, for instance spin 
or isospin/hyperfine states. E.g., the su(2) model has two internal states, which can 
be interpreted as | and |. The coordinates oi only contribute to the symmetry and 
do not enter in the functional form of the solution. 

The free parameter c acts as an attractive (repulsive) interaction strength for 
c < (c > 0) between two particle space coordinates Xi and Xj . The eigenfunctions ^ 
of Hamiltonian H and their corresponding eigenvalues have to satisfy three additional 
properties to be a solution for the quantum gas: 

(i) Symmetry, with 'J — > ivt under any 'particle exchange' (xj, Oi) <-> (xj, Oj\ i,j = 
1, . . . , TV. Here the upper (lower) sign denotes a gas of bosons (fermions). 

(ii) Continuity of $ in the interval [0, L] in each space coordinate Xi. 

(iii) Periodic boundary conditions ^(. . . x, = . . .) = V E'(. . . x, = L . , .) or hard wall 
boundary conditions , 5(. . . x = . . .) = *(. . . Xj = L . . .) = 0, i = 1, . . . , N, for 
a given box of size L. 

This model is known to be Yang-Baxter-integrable and its Bcthc Ansatz solution 
is described in detail in | Appendix A| for the given boundary conditions. Each energy 
eigenvalue is described by N' sets of complex parameters, where all these sets together 
form the coupled algebraic equations given below - the Bethe equations. For all 
models of this type the energy eigenvalue is determined solely by the quasi-momenta 
{fci, . . . , fcjv}, which are the roots of the 1st level or 1st set, with 

N 

E = J2 k l (2) 

i=l 

To determine the allowed values for {ki}, i.e., the energy eigenvalues, it is necessary to 
find a solution that satisfies the complete set of equations. Consider first the periodic 
case. These equations are 

where the formula on the left (right) is for fermions (bosons). Both fermions and 
bosons are subject to the same higher level Bethe equations, namely 

TT A i ~ A j +2 1C _ TT A i A J + 1C TT A i A j 2 lC ( .s 

11 A (/) _ A (z-i) _ i. 11 A (i) _ A (j) _ • 11 A (0 _ 1. W 

3=1 A i A j 2 1C jjti A i A j 1C 3=1 A i A j + 2 1C 

The Bethe equations for the hard wall case are 
M h - Af 5 - iic ki + Af ] - iic 



e 



2 ki L TT kj kj -f- ic kj + kj + ic fr h - kf - iic h + Aj 1} - iic 

11 fc, - kj - ic fc 4 + kj - ic 1 = 1 fe . _ A (D + i ic fc . + A W + i ic 
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where now the formula on the 1st (2nd) line is for fermions (bosons). Both bosonic 
and fcrmionic systems arc subject to the same higher level Bethe equations, 

M Af - A?-) - lie Af + Af 1 ) - |ic ii Af - Af - ic A f + Af - ic 
Yf 1 Af-Af^-jic Af+Af^-jic 
X f = [ Af -Af^ + iicAf +Af 1 ) + iic' (6) 

The above notation for the periodic and the hard wall boundary conditions 
denotes a large set of equations, in the l—th level root set {Af } there are i = 1, . . . , Mi 
roots and as many equations. To make the notation compact we have set A^ = ki and 
Mo = N as well as Mn> = 0, i.e., the product rif=a' = 1 is unity. Here the numbers 
of roots in each set is restricted by > Mi > ... > Mw-i > 0. The connection 
between the number of spin coordinates er, , . . . , er/v in the internal states 1, . . . , N' on 
the one side, and the quantum numbers Mj on the other, is given by Mjy-i = Nn>, 
Mv-2 = iViv-i + Niv, • • •, Mi = JV 2 + . . . + N N >, N = Ni + N 2 + ... + N N ,. Here 
we use a convention to label the internal state, in which the most cr^ are in, as 1, with 
N\ the number of particles in 1 . Similarly the internal state in which the second- most 
number of cr; are in is named 2 with N 2 particles in it, until we reach TV', the internal 
state with the least number of o% in it, namely Nn>. This order is no restriction on 
generality, all other values of the internal states {o~i...o~]\[} can be obtained from that 
solution via renaming the internal states accordingly. Note that Mi = is explicitly 
allowed. Thus one can obtain solutions with lower number of internal states from 
general ones with higher number of internal states. For illustration, the example of 
the reduction of spin-^ solutions to spinless solutions is shown below. Thus to obtain 
the energy eigenvalues of these models a set of N + Mi + M 2 + . . . + Mjv coupled 
equations for N+Mi+M 2 + . . .+Mjv< variables has to be solved. As already remarked, 
only the first N variables of the first level of the solution are needed to determine the 
energy eigenvalue. 

For the rest of the paper we restrict ourselves to the case N' = 2, namely two 
internal states, which can be interpreted as spin-| and spin-|. The explicit Bethe 
equations for this case reduce to four systems, hardwall bosons, periodic bosons, 
hardwall fermions and periodic fermions, which arc (from top to bottom) 



.4 



Q i2kiL 



■^L ki — Aj — |ic ki + Aj — |ic y-r ki — kj + ic fcj + kj + ic 
h — Aj + |ic ki + Aj + |ic j-J: h — kj — ic fcj + kj — ic 



M 



j-j ki - Aj - -j-j- ki - kj + ic 

1 ki Aj -(- T^ic . . kj k; ic 

j=l J 2 j^Ai J 

1 ki — Aj + |ic ki + Aj + |ic 



n 



e i2fe;L _ TT 1 3 



3=1 k i ~ A j - |ic h + Aj - |ic 
M h-A- < 1 - 



e ifc,L 



II I -Z -I c (?) 
This set (i = 1, ... , N) of N 1st level Bethe equations is connected to a second set of 
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so called spin roots {A 1; . . . , A M } obeying the M second level Bethe equations 



n J Ai — kj — iic Ai + kj — iic y-r Ai — Aj — ic A. t + Aj - 
A . _ !• . _1_ ii/> A ■ -J- I- . _L L\n ~ 11 



A l - kj + |ic Aj + % + ±ic ^ Aj - Aj + ic Aj + Aj + ic 

N » , i. M . , 

-p. Aj - k 3 - ^ic n Aj - Aj - ic 

llA.-^ + iic U A, A, • ic' U 

The 1st equation is for hard wall boundary systems, the 2nd equation is for the 
periodic models. These 2nd level Bethe equations hold for both fermion and boson 
gases, with the sign from ^ — > ±5' under particle exchange entering only in the first 
level (J7J , as can be seen from the details of the exact solution in |Appcndix A| Roots 
within one set have to be pairwise different and are usually complex numbers. Here 
M < N/2 denotes the number of one of the spin components, e.g., iV-j- = M and 
Ni = N — AI . From the above definition it can be seen that the equations for a gas 
with n internal degrees of freedom can be obtained from the solution with n' > n 
internal degrees by taking states without values in these components. For example, 
the above Bethe equations Q and JSJ| collapse for M = (blank all terms containing 
A's) to 



Y2kiL 



kj kj ~\~ ic kj ~\~ kj 



ntj — fvj -f IC fvj ; -|- h,j -f IC i2fc L i 
-J> _ -J. gi^^ _ I 
kj — kj — ic kj + kj — if 



N 



ikiL t kj — kj + ic 



— ni^f^ ^=1. (9) 

rii'i Ay-i IC 

These arc, respectively, the well known results for spinlcss bosons 8 (left) and free 
fermions (right) with hard wall boundaries (1st line) and spinless bosons £Q and free 
fermions with periodic boundary conditions (2nd line). In each equation i = 1, . . . , N . 

These equations are important, because the ground state of a Bose gas is fully 
polarized [20], i.e., described by these spinless results QJISj. The terms denoted by 
A in equations J7J and JSJl are the new terms for the hard wall case. They lead 
to the invariance of the Bethe equations under the exchange of any root kj —kj, 
A[ k ^ — > — A^ fc ' for all states of the system. From a physical point of view this means 
that the total momentum of the system is no longer conserved for hard wall boundaries, 
because reflection at the wall can change the direction but not the energy of a particle. 



3. Weak coupling expansions 

In this section we discuss the spin-i fermionic system in the weak attractive and 
weak repulsive coupling limits via expansion of the Bethe equations. For simplicity 
we concentrate on the ground state for a given polarization M, but many results and 
ideas hold as well for excitations. An explicit analytic expression for the ground state 
energy is calculated from the linearized Bethe equations. 

From numerical analysis of the Bethe equations it can be seen that the root 
distribution evolves continuously from the free fcrmions/bosons case c = 0. This 
implies that the momenta are still very close to the free momenta kj = ^r-ni, Tij = 
0, ±1, ±2, . . . (periodic) and kj — j^m, n, = 1, 2, 3, . . . (hard wall). Note that fc, = is 



Periodic and hard wall fermions and bosons 



G 



periodic 

X X 



attractive 
X X x 



~i i i i i i i r 



periodic 



X*< X0< X«X >«x >«x 



repulsive 



-871 -671 -471 -271 271 471 671 871 -871 -671 -471 -271 27C 4jt 67t 871 



hard! wall 

X X X X X 



XXX 



_I I I I I L 



attractive 



hard wall 



x«x x»x x»x x»x >»x 



_l I I I I I I L 



repulsive 



171 271 371 471 571 67t 771 871 97t 171 271 371 471 571 67t 771 87C 97t 



Figure 1. Generic root distribution in the complex plane for the ground state 
for weak attractive (c = —1, left) and weak repulsive coupling (c = 1, right). Here 
N = 14, M = 5 and L = 1, corresponding to five pairs of fermions and four 
unpaired fermions, with N\ = 9 and N-t = 5. There is a continuous evolution 
from free particles at c = 0, which are at multiples of ±%- for the periodic case 
and ~ for hard walls. Circles denote A; and crosses denote fc;. For each case the 
roots arc obtained by solving the Bethe equations Q and 151 . 



not a free (c = 0) solution in the hard wall case, due to lack of translational invariancc 
of the Hamiltonian, which has consequences for setting up integral equations for the 
ground state, cf. equations and (J2HJ- For spin-^ fermions at most two fermions 
can share one of these 'lattice' sites, thus it is natural to interpret them as a pair of 
spin-f and spin- J, . Each lattice site can be either unoccupied, occupied by one kj or 
occupied by two quasi-momenta, in which case they are interspaced by a spin root A 
(for the ground state only). In the attractive case these quasi-momenta form a complex 
pair kj = a ± i/3 with real part a almost the same as the associated real A. A typical 
ground state root distribution is shown in Figure ^ The energy k \ + k\ = a 2 — 1 
of such a pair is lower than the energy of the same two fermions in the limit c = 0, 
suggesting the notion of 'bound pair'. In the repulsive case the quasi-momenta are 
still interspaced by A, but lie on the real axis. Although there is no binding energy 
gained here, we will still use the term 'bound pair' in the following to denote a pair 
of quasi-momenta kj associated with the same site ^y-Tii or Jriz- 

The expansion of the Bethe equations follows similar work on systems containing 
bosons, fermions or mixtures of both in Batchelor el al j2SJ|3niE]. The free parameter 
in the Bethe equations is cL, but for physical clarity we will expand in small c, thinking 
of L as a fixed and finite number. Before expanding the product and the exponential 
it is necessary to ensure that i > c in /(s) = by separately treating terms 

involving quasi-momenta and spin roots from the same lattice site ^f-ni. Numerical 
analysis readily establishes the Bethe roots to be of the form 

,<■•> = ^,„, + i <«> c , {m , } _/ ± £z£zi 
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for i = 1, 



, N - 2M 



for 



1, 



Xp) 



M and 



A 



O) 



2/r 



2tt 



{ni} 



M - 1 



M — 1 



M — 1 



M- 1 



Here superscripts (p) and (u) denote paired and unpaired roots with the real constants 
/3i,7j, Si to be determined. For simplicity we consider only the case TV even and M odd, 
thus mi and rii always denote integers. In the above the ground state configuration 
for a given polarization M has all paired momenta lying symmetrically distributed 
closest to the origin with symmetrically packed unpaired momenta on an outside layer 
(recall Figure QJ. 

Assuming this functional dependency of the Bethe roots to lowest order in c it is 
possible to expand the Bethe equations consistently in powers of ^/c, leading to a set 
of linearly expanded equations to first order in c. For the periodic case this leads to 
explicit expressions 

, (u) 27T . V> 1 

k l = -r m i + c ■ 



, (r>) 27T /C 



' 2ir(mi - nj) 

M 



E 



i 



N-2M 



2-K(Ui 



J 2n(ni - rrij) 



(10) 



for the location of the Bethe roots. Similarly, for the hard wall case 

M 



IT C \ ^ 

-m, + — y ^ 



2ir 



-rii ± 



2L 

N—2M 

E 



i 



i 



c 

2^ 



M 

E 



1 



(11) 



Note the 1/ni term, which is a feature of the hard wall case. This correction has a 
significant effect and is also found in the equivalent Hubbard model expansion 
It contributes to the difference in energy for periodic and hard wall models in the 
thermodynamic limit, the so called surface energy. 

In both cases the binding strength, especially in the attractive regime, is the same 
for all pairs, which is found to be true numerically only for weak coupling ;S : . For 
stronger coupling the inner lying pairs appear to be more strongly bound (i.e., they 
have a larger imaginary part ) than the outer lying pairs, which is not visible in low 
order of c. It eventually leads to a break down of the weak coupling expansion for 
stronger interaction. Summing over the squares of the approximated roots fcj, the 
ground state energy to first order in c for the periodic case is given by |Hlj 



■MN 2 + (\M 2 - jz)N]+ —{N - M)M. (12) 

L 
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-2-1012 -2-1012 

interaction c interaction c 



Figure 2. Comparison for the ground state energy per particle between the 
analytic weak coupling results 1121 and 1131 and numerical values obtained by 
solving the full Bethe equations. Here N = 18 with M = 9 (left) and M = 3 
(right). For each case the energy is given without the free particle contribution, 
for ease of comparison. 



Similarly for the hard wall case 

E HW = (|) 2 [5 A 3 + (| - M) A 2 + (M 2 - M + i) N + M 2 ] 

+~2(A -M+l)M. (13) 

Figure [21 shows a comparative plot for the ground state energy for two different 
fillings in the weak attractive and weak repulsive regimes. The exact result (obtained 
numerically) is compared with the analytical approximations (|12[) and 1|13|) . 

4. Strong coupling expansions 

4-1- Strong repulsive limit c — ► oo 

Consider first the periodic case. As the repulsive strength of interaction is increased 
the quasi-momenta defining the ground state move away from the free values fc; = 
^rijjTii = 0, ±1,±2, ... towards the values fcj = Jrij,nj = ±1,±3, ... for tightly 
bound fermions. For the hard wall case the quasi-momenta return to the free values 
ki = j-rii, rii = 1,2,.... The spin roots A wander out to infinity linearly with c, except 
for at most one A = in the periodic case [TJ\ . For the hard wall case the behavior is 
similar. An illustration of the root behavior for the different cases is given in Figure|21 
For the ground state all roots remain real. The energy is given in terms of a set of M 
parameters {7^}, appearing as the leading coefficients in the spin roots Aj. 

In this limit, the roots are given in integer powers of c, in agreement with 
numerical checks. For the periodic repulsive case the ground state roots are given 
by (see also chapter 2.27 of Ref. (17) 1 

k i = Ym + kc~ 1 i = l,...,N {n q } = {±1,±3, . ..,±(N- 1)} 
Aj = 7iC i = 1, . . . , M (with one A = for M odd) 
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Figure 3. The ground state quasi-momcnta distribution for periodic and hard 
wall fermions in the strong coupling regime c — » ±oo for different polarization M . 
(top left) Strongly attractive periodic case. Complex conjugate pairs of quasi- 
momenta k axe approximately interspaced by spin roots A with k = A ± i^c. 
(bottom left) Strongly attractive hard wall case. Similar to the attractive periodic 
case, (top right) Strongly repulsive periodic case. Here the spin roots A wander 
out proportionally with c towards ±co but remain real 1171 . The quasi-momcnta 
ki break up the 'bound pairs' and tend to their tightly bound fermion values 
ki = j^rii,ni = ±1,±3, .... (bottom right) Strongly repulsive hard wall case. 
Similar to the repulsive periodic case. The spin roots A wander to +oo, but the 
quasi-momenta tend to the tightly bound fermion values k{ = j^ni,rii = 1,2,... 
breaking up 'bound pairs' (see text). In each figure the quasi-momenta k and 
the spin roots A are denoted by crosses and circles where applicable. The 
numerical values are N = 8,M — 3,c = 500, L = 1 (repulsive cases) and 
N = 14, M = 5, c = —200,L = 1 (attractive cases). The dotted lines arc guides 
to the eye. 



where we show only orders which are relevant for the energy to order 



The 



real constants Si and Aj are to be determined. Expanding the second level Bcthe 
equations up to order c _1 the parameters 7, are found to satisfy the same parameter- 
free equations 

N M 



H + 2 1 



n 



7i 

H ~ li + - 



for i = l. 



,M 



(14) 



as for the isotropic Heisenbcrg spin chain. From the expansion of the first level Bcthe 
equation up to order c _1 the roots are found to be (for i = 1, . . . , M) 



7T 

— Uj 

L 



1 - — n 



cL 



Note that the constant fl^.j 



(15) 

depends only on the set {7i}i=i,...,M 



J=l 7?+l/4 

and neither on the interaction c nor on the integers n^, making numerical evaluation 
significantly easier than for the initial Bethc equations. 
The energy expression up to order c _1 follows as 



1 n 

cL 



0{c 



(16) 
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Figure 4. Comparison of the ground state energy in the strong repulsive 
interaction limit c — » oo for (left) periodic boundaries and (right) hard walls. The 
circles denote the values obtained from the leading terms evaluated in equations 
1161 and 1191 , The full curves are the results obtained from the exact numerical 
solution of the Bcthe equations. The inset shows confirmation of the order c 
coefficient. Dashed lines denote numerical values obtained from (E num — S aaymp t)-c 
while full lines are the c~ 1 coefficients from equations 1161 and 1191 . The lines 
converge to the same value for large c. The four cases shown are as in the main 
panel. The same numerical values are used as in Figure |2] namely N = 18 and 
M = 3,9 with L = 1. 



where we have taken the integers rij for the ground state configuration. Comparisons 
between this expansion and the exact numerical solution for finite systems are given 
in Figure |U Formulas for the hard wall case are derived in a similar fashion. Here 



hi — ^ rii 
A, = 7i c 



Sic 



M 



{m} = {i,2,...,N} 



in the strong coupling limit c — ► oo. Expansion of the second level Bethe equations 
results in the spin root coefficients ji satisfying 

2N M 



li 



li 



n 



li - 1] + i li + 7j + 1 
1% - lj - i li + lj - i 



(17) 



for i = 1, . . . , M. These are the well known Bethe equations for the open Heisenberg 
chain Expansion of the first level Bethe equations to order c _1 gives the 

quasi- momenta {ki} depending only on the parameters {7^} via 

1 



Here the constant Cl^j 



1 



cL 



E 



M 
3=1 



7? + 1/4 



is obtained from the solution {7i}j=i, 



equations l|17|) . The leading terms in the ground state energy follow as 

^2 



Ehw = 



6L- 



(2N 3 + 3N 2 + N) 



i-i^yv 

c L 



M 



0{c 



(18) 

,M Of 

(19) 



A comparison for the ground state energy between this expansion and the exact 
numerical result is shown in Figure 0] for a finite system with different polarizations. 




Figure 5. Ground state energy per particle for repulsive periodic and hard wall 
fcrmions versus interaction strength c. The dotted lines denote the values for 
free c = and tight binding c — » oo fcrmions. Full polarization corresponds to 
non-interacting fcrmions and is thus only a horizontal line (not shown) . For given 
polarization the periodic system always has a lower energy than the corresponding 
hard wall system. Shown are numerical values for N = 18 and M = 3, 6, 9 with 
L = 1. 



The results are in excellent agreement with the numerical curve as soon as the 
assumption A ~ c holds, which happens for certain system sizes N and M at large 
interaction strength c > 1000. Note that here we have plotted against the interaction 
strength c, one can also do similar plots against the rescaled interaction 7. 

Figure [5] shows a double logarithmic plot of the ground state for a generic 
finite system with different polarization, obtained numerically from the original Bcthe 
equations. The dotted lines denote the known asymptotic values for free (c — ► 0) and 
tightly bound fcrmions (c — > 00), which agree with the leading order of the expansions 
derived here. Note that for a finite size system the hardwall case has a higher energy 
than the periodic case. 

4-2. Strong attractive limit c — > —00 

In this section we again for convenience consider the case N even with M odd, the 
generalization from which should be clear. The root configuration for the ground 
state consists of 2M complex conjugated quasi-momenta k, each pair interspaced by 
a real spin root A. The remaining N — 2M unpaired momenta k lie on the real axis. 
The occurrence of complex roots makes the analysis more complicated. Following 
Takahashi we assume fcj = Aj ± ^ic for the 2M paired quasi-momenta. This is a 
good approximation for the thermodynamic limit |17j . We assume here that it is also 
a good approximation for large attractive interaction strength c — > — 00. Using this 
form the Bcthc equations can be rewritten (see chapter 2.2.8 of Ref. ^3) m terms of 
the root sets {ki}i=2M+i,...,N and {Ai}i = i t ...^M, i.e., consisting only of the remaining 
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(real) quasi- momenta and the (real) spin roots Aj, as 



A I 



e ik iL = 



■j-j- ki - Aj + 5 ic 



N 



A I 



A2A,L 



A, 



n A t kj + |lL yj 

, A,- - «„■ - hie -14 A 



A, - Aj + 



2 1C 



A,- 



12 



(20) 



All results in this subsection refer to the analysis and numerical evaluation of these 
modified Bethc equations. Asymptotic and numerical analysis suggests that the roots 
have the form (n, , m, G Z) 



«i = m-f H — , 
L c 



A; 



7T 



7i 
c 



2M+1,. 
l,...,iV, 



1 A r , 



{±1,±3,. 
M- 1 



,±(iV- 2M- 1)} 
M-l' 



Expansion in the small parameter c 1 results in 



h = rrii— 1 - 



Aj = rii- ( 1 - 



4M 

~Lc~ 

2N- 3M 
Lc 



(21) 



This root distribution gives the energy 



A I 



E 



PBC 



E K 2 



2 
L 2 



N 



2 1 



E * 2 

i=2M+l 

4iV-6M\ M 3 - ill 



Lc 



12 



8M 
~Lc~ 



| [iV 3 - 6MN 2 + (12M 2 -\)N + 2M- 8Af 3 ] | 



0(c 



(22) 



The derivation for the hard wall case follows the same steps, with the 
corresponding modified Bethc equations 



A4AiL 



M 

n 



ki — Aj + iic ki + Aj 



_ ± ki — Aj — ^ic ki + Aj — iic' 



2M + 1, 



n 

j=2M+l 



Aj — fej + iic Aj + kj + iic 



A, 



ic A, 



A; 



Aj + ic A* + Aj + ic 



,-ic " Aj 



A, 



ic A l + Aj - ic 



(23) 



for i — 1, . . . , M. Expanding these equations to order c 1 leads to 



A, 



L 

L 



—M 

Lc 

2N - 3M - 
Lc 



2M + 1, 



,AT m j = {l > 2,...,JV-2M} 



1,. 



,M 



M 

Y 



Using these expressions for the roots gives the expansion for the energy 



Ehw = - \Mc 2 + — 
8M\ 



AN - 6M - 2 



Lc 



M(M + 1)(2M + 1) 



UN - 2M)(N - 2M + 1)(2N - AM + 1) 



+ 0(c- 2 ). (24) 
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attractive interaction Id attractive interaction Id 



Figure 6. Ground state energy in the strong attractive region. The curves show 
the comparison between the numerical solution of the modified Bethe equations 
121 H and 1231 [solid lines] and the expansions and @U [dashed lines]. For 

clearer comparison both the binding energy E^ ni ^ ~ c 2 of the paired quasi- 
momcnta and the O(c ) energy term is subtracted. The values used are TV = 18 
with M = 3 (upper curves) and M = 9 (lower curves) with L = 1. 



Figure [S] shows a comparison between the expansions l)22|) and 1|24|) and the exact 
ground state energy for the strong attractive case for periodic and hard wall boundary 
conditions. Note that both the expansion and the exact numerical solution refer to 
the modified Bethe equations with paired quasi- moment a of the form k = A ± iic. 

5. Integral equations for the ground state energy 

It is convenient to introduce the standard logarithmic form (see, e.g., Ref. |TJ|) of 
the Bethe equations JJJ) and JHJ. In terms of the function <j>{x) = arctan(x/c), the 
equations for periodic fermions are 

M 

fc i L = 27rI < -2 5^0(2(k-A J )) 

N M 

2 <f> (2{Ai - k )) = 2-n-Jj + 2 0(A* - A 3 ) (25) 
j'=i 3=i 
and for hard wall fermions 

M 

2kiL = 2itl t - 2^ [<t>{2{h - Ajj) + <j> (2(fci + A,))] 

N 

2 ( 2 ( A *-^)) + 0(2(A 4 + *,•))] 
j'=i 

A/ 

= 27r Jj + 2 ^ ^(Ai - Aj) + cp{A t + A,)] . (26) 

The constants I±, . . . , In and Ji, . . . , Jm act as quantum numbers, they are either all 
integer or half-odd integer within one set, depending on the boundary conditions and 
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if N and M are even or odd ^7]. Generally the solution for a given allowed set of 
quantum numbers is unique [T^]. The ground state for the fermionic system is given 
by the 'half- filling' M = N/2 [TJ\, while in comparison the ground state of bosonic 
systems of this type is totally polarized (M = 0) [2U|. As a result, the ground state 
analysis of bosonic systems can be done with the one-level Bethe equations for spinless 
bosons (given above), as long as no special polarization is required. 

We only consider here the special case of fermionic systems with N even and M 
odd. The quantum numbers I n and J n for the ground state configuration are all one 
unit away from each other and either integers or half-odd integers, as they result from 
counting phases. For excited states the quantum numbers become more spaced. For 
the ground state in the attractive case a more convenient, approximative form (see 
eqns 2.115 and 2.116 in Ref. [2j as well as (|20|) . (|23|1 in the previous section) of the 
Bethe equations in which all roots are real can be used. Analogous equations hold 
for the hard wall case. From these equations the behavior in the thermodynamic 
limit can be obtained via the standard integral equation procedure [5] El 1171 1151 IT3*| . 
In the limit N, L — > oo with constant particle density N/L it is easiest to consider 
only the densities p(k) and cr(A) of the root sets {ki} and {A^}, which are obtained 
from knowledge of the distribution of the quantum numbers I n and J n in the various 
regimes and for the different boundary conditions. 

In the repulsive periodic case the ground state roots are all real and their density 
is determined by the coupled equations 

1 2c f B cr(A)dA 



2tt tt J_ b c 2 + 4(ft - A) 2 

c f B <r(A')dA' 2c [Q p(k)dk 
1 ' ~ tt J_ B c 2 + (A — A') 2 7r J_ Q c 2 + 4(fc - A) 2 {Z ° 

The density n and the polarization n± determine the constants B and Q in the 
equations below, making it possible to calculate the energy per particle e, with 

rQ pB 
k 2 p{k)Ak 7 n = / p(k)dk, n i = / (r(A)dA. (28) 



J-Q J-B 

For the attractive case, the possibility of bound pairs changes the form of these 
equations to 

° B a(A)dA 



1 


2c 


2^ 




1 

- -f 




TT 


71" J-B 



3 c 2 +4(/c-A) 2 

cr(A')dA' 2c f Q p(k)dk 



(29) 



c? + (A - A') 2 tt J_ Q c 2 + 4(fc - A) 2 

with the auxiliary conditions 

fB fQ fQ 

n = 2 cr(A)dA+ / p(k)dk, n l - n T = / p(k)dk 

J-B J-Q J-Q 

r 2k 2 - -c 2 ) cr(A)dA + / k 2 p(k)dk. (30) 



—B 



2 



These equations have been extensively examined in the periodic case (see, e.g., 
|37l 1ST] and references therein) . To obtain the integral equations for the hard wall 
case it is necessary to proceed with more caution - the lack of translational invariancc 
forbids the roots to take the value zero, even though it leads to a solution of the 
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Bcthc equations. It is convenient to symmetrize the exponential Bethe equations 
to include sets of zeros of the form {±£^,0} and {±Aj,0} ^H2HI- This makes the 
equations somewhat resemble a special system of 2AT+1 particles (see also Ref. ^3] f° r 
a remark on the equivalence of spinless bosons with hard walls and periodic boundary 
conditions). The procedure then becomes analogous to the periodic case (e.g., see 
chapters 2.2.6-2.2.8 in Ref. [Uj and J2H1) - It has, e.g., been applied to the XXZ 
chain j^HI and the open end Hubbard model |28l I38j . Note that the auxiliary integrals 
contain additional factors and terms accounting for the zero value. The resulting 
integral equations for repulsive interactions and hard walls are 

1 2c f B cr(A)dA 1 2c 1 

p{k) = - 



TT 



l_ B c 2 + 4(fc - A) 2 L TT c 2 + 4fc 2 
c f B cr(A')dA' 2c [Q p(k)dk 1 c 1 

1 ' TT J_ B C 2 + (A - A') 2 7T J_ Q c 2 + 4(fc - A) 2 L TT c 2 + A 2 1 ' 

with the restricting conditions 

I rQ i i rB 1 

«=o/ P(k)dk-—, n t = - cr(A) dA — 



2 J q 2 J q *2iL 

i 

2 



k 2 p(k)dk. (32) 



This agrees with the result for the general su(N) fermionic case when the 

boundary potentials in their model become infinitely strong and the extra terms 
vanish. The same calculation for attractive interactions and hard walls leads to 

1 2 f B ccr(A)dA 2 c 



p(k) = 



TT 



_ B c? + 4(k - A) 2 ttL c 2 + 4fc 2 



m = 1 + I [ B c<7(AQdA^ 2 /■« cp(fc)dfc 
CTl ; - + ~ l_ B c 2 + (A-A') 2 W-n c 2 +4(A-fc) 2 



7T 

1 c 2 



ttL c 2 + A 2 Tri c 2 + 4A 2 
with the restricting conditions 

i /■« ,, 3 i /■« i 



(33) 



2 j ^p(k)dk+ J ^a(A)dA-— , n; -n T = -y p(fc)dfc 2£ 



Q 

fc 2 i o(fc)dfc+ / (A 2 -c 2 /4)ct(A)c1A. (34) 

B 



In these equation terms proportional to are 'surface energy' terms, which appear 
for the hard wall case compared to the periodic case. The L appearing after taking the 
thermodynamic limit is understand to be large but not infinite. The only effect of the 
different boundaries appears for this still finite L. Apart from these terms the integral 
equations for the attractive and repulsive case are the same as those for the periodic 
case. The above hard wall integral equations can be used to obtain expressions for 
the surface energy and then compared to the results obtained by taking a continuum 
limit in the Hubbard model \27\ . 

Note that in the derivation of the ground state integral equations in the attractive 
regime bound pairs of approximate form k a = A a ± ^ic were assumed (see also eqn. 
2.114 in Ref. ^7j). This is no contradiction to the dependency of the bound momenta 
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on the coupling strength going as y/c in Section The equations are valid in 

different regimes, for weak coupling we examine a finite size (i.e., L and N finite) with 
vanishingly small coupling constant c. In contrast, in this section we consider a system 
at fixed finite c and take the thermodynamic limit N, L —> oo with N/L = const. These 
limits do not commute, see the discussion in Rcf. [2SJ 

6. Concluding remarks 

In this paper we have extended the Bethe Ansatz solution for one-dimensional bosons 
and fermions with S- interaction and arbitrary internal degrees of freedom to the case of 
hard wall boundary conditions. The general solution Q and © relies on the known 
Bethe Ansatz solution of the open su(N) chain. The special case of the two-state 
fcrmionic model has been analyzed in detail, both for periodic and hard wall boundary 
conditions. Expansions were obtained for the ground state energy in the weak and 
strong coupling limits in both the attractive and repulsive regimes. Integral equations 
were also given for the ground state. We expect that given the recent experimental 
advances on one-dimensional quantum gases the exact solution to the family of models 
considered here will enjoy further interest, particularly for the investigation of finite 
systems with a boundary. 
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Appendix A. Coordinate Bethe Ansatz for ^-interaction quantum gases 

In this Appendix we give the coordinate Bethe Ansatz solution for the <5-interaction 
gas with hard wall boundary conditions. For case of comparison we also discuss the 
more standard case of periodic boundary conditions. Our aim is to give a self-contained 
treatment. For further background the reader is referred to the books [1511171 171 11311141 
and where appropriate, to the original literature. Like for the periodic case, the 
problem is solved by making contact with results for the diagonalisation of the su(N) 
chain, the solution of which is discussed in |Appcndix B| 

Unlike the spin chains, the Hamiltonian does not have a mechanism for 
spin flipping. However, the model is block diagonal for different spin arguments 
{Ni, . . . , Njyi}, where JVj is the number of coordinates <7j having value i. For example, 
for spin-i particles the number of spin f and spin J, particles is conserved. We want 
(anti-)symmetry under coordinate exchange, as we are considering indistinguishable 
particles. The 2-particle ^-interaction scatters elastically, i.e., at most momenta can be 
exchanged, but not their value (however, scattering at a reflecting wall might change 
fcj — * — ki, as discussed below). The contact interaction 5(xi — Xj) only plays a role 
'near' the hyperplanes, where two coordinates are the same Xi — > Xj. In quantum 
mechanics books j^U El it is shown by integrating over the region boundary, that 
the wave function is continuous but has a finite jump in the first derivative. Thus 
generally these problems can be seen as simple free particle problems away from the 
hyperplanes Xi — > Xj, with plane wave solutions. These regions with free solutions are 
connected by appropriate boundary conditions, accommodating the jump in the 1st 
derivative from the contact interaction. This consideration leads to the Ansatz wave 
functions 

. . .,x N ,ax, . . . ,<t n ) 

= J2 A <nr:°APi, Pn\Qi, Q N )e i ( kp ^+-+ k ™ x ^ (A.l) 
p 

^{xi, . . .,x N ,ai, . . . ,tjjv) 

= E E ^....^(epi, . • ■ , t PN \Qi, Qiv)e i(epifepia:Q1+ - +epJvfepJVXQJv) (A.2) 

P cpi=± 

for arguments Xi ^ Xj Vi ^ j with i,j = 1,...,N. These wavefunctions are for 
periodic and hard wall boundary conditions, respectively. 

In (|A.lfl and (|A.2fl the sums run over all AH permutations P of {1,2, . . . , JV}, 
and over all N signs ti = ±. The wavefunctions are piece- wise defined functions - 
for each given argument (x\, . . . , Xn) there is exactly one permutation Q, so that 
< xqi < XQ2 < . . . < xqn < L. For (jA.lj) there are N complex, pair-wise numbers 
{ki} and Nix Nix N' N complex numbers A ai) ___ j(TN (P\Q) labelled by the permutations 
P, Q of the combinations and signs. Similarly for i|A.2[> there are JY! x JV! x N' N x 2 N 
complex numbers A ai o- n (^pi---^pn\Q)- 

We turn now to the key properties of the wavefunction. 
(Anti-) Symmetry. Applying the (anti-)symmetry conditions ^ — ► ^ (bosons) and 
^ — > — ^ for (fermions) under particle exchange [xi,Ui) <-> (xj,o~j) to the wave 
functions (jA.ljl and (| A. 2|) leads to (For brevity these formula are for hard walls, the 
periodic case is recovered by setting epi — > Pj.) 

A ai aj (e Pl , . . . ,ep N \ Q a Qb ) = ±A aj CTl (e Pl , . . . , e Pjv | Q b Q a ) (A. 3) 
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where we have introduced an obvious abbreviation in the notation. Here Q a and Q b 
denote the positions of Xi and Xj for a given argument of and we assume Xj < Xj 
without loss of generality. In the same way define i — P a and j = P b . Introduce the 

operator T» via (T«)^;"^ = ±6 <n,^ 6 ^,<r[ Xl r ?ij ^,<> i-e., T« = P« for bosons 
and T y = — T 5 *- 7 for fermions, in terms of the permutation operator V %3 . The (anti- 
) symmetry condition can then be rewritten, using the summation convention for 
internal matrix indices from now on, as 

( t<j ')^:::^ 1 ...^(^ > ...,6p,iq qo 

= A ai ... (JN {ep 1 ,. ..,ep N \ Q b Q a ). (A.4) 

Continuity. At the intersection of two regions the piece-wise definition of the wave 
function ^(xx, xjy, a\, crjy) *%<*) = ^(xi, ■ ■ ■ , %n, a%, . . . , crjv) *<>*i , requires 
that the coefficients satisfy 

A ai ...a N { £P a £P b I Qa Qb ) + A tTl ... <TW ( e P 6 e P„ I Qa Qb ) 

= A . 1 ... . JV ( ep a e Pb \ Q b Qa) +A a . 1 ... a . N ( E Pb e Pa I Q b Q a ). (A.5) 

5 -interaction. In standard QM textbooks (e.g., |34l I14j ). it is shown by integration 
that in local center-of-mass coordinates X, y for the two coordinates in the argument 
Xi = Xg a — > Xj = XQb the (5-interaction leads to 



dy 



y=0+ d V 



= c^\ y=0 . (A.6) 

v=o- 

Applying this last result to the wavefunction results in a relation involving four 
A coefficients with different permutations P and Q, one of which can be eliminated 
by using (|A.5() . In this way 

i(k Pb - k Pa ) [A az {j. ( P a P b I Q a Q b ) - A a . Vj ( P b P a I Qb Qa )] 

= C [A ai 0j (P a Pb\Qa Qb)+A ai CT . ( P b P a \ Qa Qb )] . (A.7) 

This equation, connecting A coefficients with different permutations P and Q but 
identical spin indices o~i, can be transformed, using the operator T from (|A.4f) . into an 
equation connecting A with a different permuation Q, but the same permutation P 
and spin indices <jj by 'hiding' the spin coordinate permutation inside the definition 
of T %J . The resulting central equation was found by Yang [5]. In our notation it takes 
the form of a matrix equation 



A ( j 1 ... < j n ( P a Pb I Qa Qb ) 



i(k Pb -k Pa )T l i +cl& 
i(kp b -k Pa ) - c 



Aa[...a' N ( Pb Pa \QaQb) 

ai...a N 



in the spin indices o~\ . . . un . Here the highlighted term defines the operator 
Y l3 (kp b — kp a ). The formula for hard walls is obtained by replacing fc^ — > ejfcj and 
also Pi — > epj in the coefficients. For future reference we introduce the operators 

\X i i(u)V' 1 "' a ' N := rr y '(u)li/''^ \T i Uu)Y"-' T " , which connect coefficients with 
permutations P and Q differing by the same transposition only 

M-'J P - P b\QaQb)= [X ij (k Pb - k Pa )}°};^l A a ,... a , N ( P b P a \QbQa) 

So far restrictions have only been imposed on the coefficients A, while the N complex 
ki are arbitrary (the 'scattering problem'). Imposing boundary conditions restricts 
the choice of these values to the so called Bethe equations, which we now turn to. 
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Appendix A.l. Periodic Boundary Conditions 

Without losing generality the periodicity condition can be written as ^l^^o = ^Ui=L) 
which when applied to the ansatz wave function IjA.lfl leads to the condition 

A ai ...a N {Pl, P2, ■ ■ ■ ,Pn\Qi, ■ ■ ■ iQn) 

= e ik ^ L A^,„ aN (P 2 , ...,P N , P x \Qi, Q N ). (A.8) 

With the help of repeated application of the X operators defined above this can be 
rewritten as an equation connecting coefficients A(P\Q) with identical permutations 
P and Q through a series of matrix multiplications. The way of commuting Qi and 
Pi through needs to be path- independent to be consistent. This is ensured by the 
Yang-Baxter equation, which is satisfied by the X operators. This can be seen by the 
identification of the R and X operators in 1A.10|) below and then using 1|B.2(1 . 
Applying successive multiple X operators leads to the set of N equations 

A ai ,.. aN (l...i-l,i+l...N,i\l...i-l,i+l...N, i) c ik ' L 

= X u (h - ki)f. . . X 4 ~ l4 (fc 4 _i - hf 3 lX u+1 (k t+1 



X N \k 



N 



A a 



(l...i-l,i+l...N,i\l...i-l,i+l...N, i) 



which can be simultaneously diagonalized. Here we have introduced the abbreviated 
notation a = (<7 1 ...erjv). Note also that P-2 = 1, ...,Pjv = N an d wc choose 
P\ = i and adapt the indices accordingly. This means it is possible find linear 
combinations of Ag(P\Q) considered as vectors in the indices ai ...ctjv- From now 
on we suppress the permutations, as it is always the same argument (P\Q) in the 
remaining equations. Call these linear combinations say, MS A-,, so that the right 
hand side product of X operators can be replaced with its eigenvalue 9{ki\k\ . . . /cat) 



via e 



kiL 



MA = X ... XMA = 6{. . .)MA, resulting in the 1st level Bethe equations 

e ifc - L = 6(k l \k 1 ...k N ) Vi = l,...,N. (A.9) 

The eigenvalue 9 can be obtained on rewriting the product of X operators in 
terms of R matrices via the identification 



X b {u) = - — ^Riiu/c) X f (u) = R(-i u/c). 



(A.10) 



u + 1 c 

Note that the identification for bosonic and fermionic systems is different. The R- 
matrix is defined in (|B.1|) . The X operators contain a signed permutation operator 
for either type of particle, resulting from the condition of (anti-)symmctry under 
particle exchange (we suppress matrix indices). Thus 



X U (k! - h) . . . X 1 - 11 ^! - h)X 



li+1 



n' kj - kj + ic 
ki - k, - ic 



tr f 



Ri 



(ki+i - h) . . . X l (kN-h) 
(u-ki)) ... R N0 ( --(u-k N ) 



tr 



i? 



10 



■(u - ki) 



R 



Ni) 



few) 



(A.ll) 



u— ki 



Here we have added a trace operation over an additional space of the same type as 
the other N spaces. The second and third lines in this equation are for bosonic and 
fermionic X, respectively. 

Because the R matrix reduces to the permutation operator at Rio(i(ki — ki)/c) = 
P t0 , the trace can be evaluated and this product reduces to a product of the X 
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operators. The operator in the square bracket is the monodromy matrix of the periodic 
su(N) chain. Its eigenvalue was originally obtained by Sutherland [7]. We give the 
resulting eigenvalue and also sketch the technically more interesting diagonalisation 
of the open su(N) chain, of relevance to the hard wall ^-interaction quantum gases, 
in |Appcndix B| 

Using the eigenvalue A(u) from equation (|B.3|) gives the 1st level Bcthe equations 
l|A.9|) in the form 

T -p r A I u = — h {q,} = {--kj} (bosons) 



e ik,L 



e ife * L = A ( u = +-ki 



{ij} = ) (fermions) 



e ifc,L 



c 

for i = 1, . . . , N. Using a shift in the roots gives the standard form 

fl f ~^T? iC (bosons) (A.12) 

Ki Kj 1C Ki A-i -f- ttIC 

M Z — A _i_ 1 ■ 

n r C ( f — ) (A.i3) 

The momenta fci, . . . , fcjv are constricted by A/ additional Bcthe equations (here 
for su(2), the general su(N) case is given in equation Q) 

The energy eigenvalue associated with the Bethe wavefunction is given in terms of the 
1st level Bethe roots, regardless of the number of nested levels, by 

N 
i=l 

because the initial ansatz was just a combination of plane waves with appropriate 
boundary conditions. To obtain the general solution for N' internal states insert the 
general eigenvalue given in i|B.3|l into the right hand side of (|A.9|) and use the N' — 1 
levels of Bcthe equations in l|B.4() . 

Appendix A. 2. Hard Wall Boundary Conditions 

The hard wall boundary condition ^I^—q = ^\ Xi= L = for i = 1,...,N, translates 
into the equivalent of l|A.8[) when applied to the wave function l|A.2|) . namely 



Ari...cr N (ep 1 =+, ep 2I ■ ■ ■ > e P N \Qii Qii ■ ■ ■ , Qn) 



— -e 1 Pl A (Tl ___ crN (ep 1 = -, ep 2 , . . . , ep N \Qi, Q2, ■ ■ ■ , Qt 



c ik ^ L A ai ... (JN (ep 2 , . ■.,ep N ,e Pl =+\Q 2 ,...,Q N ,Qi) 

= -A ai ... aN (ep 2 , . . . , ep N e Pl = — \Q 2 , ■ ■ ■ , Qn, Qi)- 

The strategy is similar to that for the periodic case described in the previous section. 
With the help of the X operators both equations can be reduced to contain the same 
permutations P and Q as arguments. Then one of the A coefficients, say the one 
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containing ep 1 
periodic case: 



is eliminated, leading to a single matrix equation S as in the 



i2fc ; L 



,e N \i,l,2 N) 



[X u (eiki - h) . ..X m (e N k N - ki)X m (-[e N k N + h}) 



(a 1 ...a N ) 



0» 

r V 6 



-,ex,...,ejv|i, 1,2, ...,N). 



This reduction is again consistent, i.e., independent of the path of transpositions, 
due to the Yang-Baxter equations. This expression can again considered to be an 
eigenvalue equation in which the product of X operators on the right hand side has 
the eigenvalue 8 HW (ki\ki . . . fcjy), leading to 

e i2kiL = e HW (k i \k l ...k N ) (A.15) 

Using the same idea as in the periodic case, this eigenvalue can be obtained in terms of 
the eigenvalue of the open su(N) chain, which is derived in |Appcndix B| at a special 
value of the argument. Taking the trace in the auxiliary space has to be done more 
carefully, in order to account for the term X iQ (— 2ki), which does not leave a prcfactor, 
because tro Rio(—i2ki/c) = 1. We now set e n = + without loss of generality. The next 
step is analogue to the periodic case, using the identification from l|A.10|) to obtain 

X u (h - h) . .. - h)X l+u {k l+1 - h) . ..X m (k N -ki) . . . 
■ ■ ■ X m (-[k N + ki})... X i+u (-[k i+1 + k^X^i-lh-x +ki))... X^i-lh 

N 



n 



tr Xb I u 



c— i(kj — ki) c—i(kj + ki) 
c+i(kj — ki) c+i(kj + ki) 



9j 



tr T I u = --h 



(A.16) 



The last two lines in this equation are for bosons and fermions, respectively. Using 
the eigenvalue A HW (it) of troXb( u ) calculated in (|B.21|) the 1st level Bethe equation 
can be written as 

M 



S2kiL 



j-j ki Aj 

3=1 

JV 

n 



gic k l + Aj + |ic 



Aj - Uc ki + Aj 



ic 



kj+ic ki + kj +ic -j-j- ki — Aj — |ic ki 



Ai 



ric 



hjj hjj \c 



kj — ic 



A, 



2 IC /l i 



Ai 



(fermions) 



(bosons) (A. 17) 



j.£i -j j j = i 

Again the 1st level Bethe quasi-momcnta {ki, . . . , kj^} are governed by a second 
set of parameters, the spin roots Ai, . . . , Am, which obey the Bethe equations l|B.22|) . 
The most general su(N) case is given in JSJ, for su(2) they reduce to 



AI 

n 



Ai - Ai + ic Aa + Ai 



N A 
IC T-r A,; 



Ai - A, - ic A, 



A, 



n 



kj 



Ai 



ic A 4 



ic A, 



^3 _ -, 



(A.16 



valid for both bosons and fermions. These equations are invariant under a simul- 
taneous change of sign of all parameters within cither set, {ki} or {A^}. 

S Again we set Pi = i, P2 = 1, . . . , Pjv = N and similarly for Q. Note that this time the term with 
X" is omitted twice. 
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Appendix B. Algebraic Bethe Ansatz for the su(N) chain 

In this Appendix we give the necessary details for the diagonalization of the su(N) 
chain, first for periodic and then for hard wall boundary conditions. Our approach 
here is via the Algebraic Bethe Ansatz (see, e.g., Ref. \Rj\), although we could have 
equally used the analytic Bethe Ansatz [391 140| . 

Appendix B.l. Periodic boundary conditions 



The diagonalization of the product of i?-operators in equations (|A.11(I was first done in 
the context of the periodic su(N) chain |7j. In the Algebraic Bethe Ansatz formalism 
conventions and prefactors differ, often making comparison between different papers 
difficult. Here we follow the normalization notation given by the rational i?-matrix 

Rij(u)= uId * j + Pij & RZUu) = -^-S atl Sp, s + -^—S a>s 5 f3 , 1 (B.l) 
1 + u M 1 + u 1 + u 

The standard notation for its entries are a(u) = 1, b(u) = u/(l + u) and c(u) = 

1/(1 + u). R acts on the two A'-dimensional vector spaces Vi and Vj, V is the 

permutation operator. Spaces with indices 1, . . . , N are denoted as quantum spaces, 

with index the auxiliary space. All relations involving i?-matrices are simple matrix 

multiplications with summations over internal indices. Proofs are most simply done in 

a graphical notation |3(>j . here we only give the final key equations here. The above R- 

matrix is normalized, with inverse (on both spaces) Rij(—u)Rij(u) = Idy. It satisfies 

the Yang-Baxter relation (YBE) 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 {u)R 12 {u - v). (B.2) 

These equations arc sufficient to prove all other relations in this section. 

For the periodic case the monodromy matrix has N arbitrary inhomogeneities 
{gi, . . . , q N } and is given by 

T (u) = Rw(u — qi)... Rno(u - qn)- 

For these conventions the eigenvalue A (it) of the transfer matrix, the 0-trace of the 
monodromy matrix T Q (u), is given by 

N Ah _ (1) _ 

A(n\{qi}) = Y[ a ( u — Qj) Y\_ (T) ^ vamsmn g terms. (B.3) 



j — j - 
=i 

Here vanishing terms denotes the remaining part of the periodic su(N) chain eigen- 
value, which vanishes due to the special argument u = qj under consideration for 
the quantum gas model. The product over the N roots is equal to unity due to the 
choice of i?-matrix normalization. In the complete eigenvalue expression all functions 
containing terms of the form IljLi b(u — qi) vanish and only the term stemming from 
the action of Tl(u) in the monodromy matrix remains (for the full expression see 
Refs. [3 El)- Generally it is most convenient to obtain the Bethe equations 

N W 11 A/ i ,,W „W 4- 1 A/ 2 ,,W ,,(2) I 

TT H V + 1 _ TT Z M J TT fc ~ ^ — 

11 (1) 11 (1) (1) , 11 (1) (2) 

i=l t'r - 1] , = 1 /',, -{*) - 1 j=l /'/. - ft) 

W-i „(0 (/— 1) . M l (0 „(/) , i Af,+i (I) 1 

TT Pfc H'j -r 1 _ tt Pfc Mj + l "i r Mfc 1 

11 (i) (i-i) ~ 11 (0 (0 -, 11 



(0 (i-i) 11 (i) (i) i 11 (i) (i+i) 
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n 



(JV'- 



1) „(JV'-2) 



(JV'- 



/' 



(N'-2) 



_ (N>- 



n 



(JV'- 



.(JV' 



- 1 



(B.4) 



from the full expression of the eigenvalue as a pole-free condition than finding the 
Bethe equations from the requirement that 'unwanted terms' in the Algebraic Bcthe 
Ansatz cancel. Here k = 1, . . . , Mi in each line, with I = 2, . . . , N' — 2 in line two. 
The standard form of the Bethe equations for the quantum gas model has a more 



symmetric form obtained by a shift in the parameters /xV"' — ► pAj 'c ' — m 
which the upper (lower) sign is for bosons (fermions) (cf. equation Details on 

performing the calculations above can be found in Refs. among others. 



Appendix B.2. Hard Wall boundary conditions 

In this section the relevant diagonalization of the double-row transfer matrix was first 
done for the su(2) case [35] in the context of the open XX Z spin chain. A pedagogical 
derivation of the general su(N) case via the Algebraic Bcthe Ansatz has also been given 
|36j . The quantum gas model requires only the special limit of the rational i?-matrix, 
with all reflection matrices being simply scalar multiples of the identity matrix. We 
will state the results and key steps adapted from Ref. j2U to the special case at hand. 

First define the monodromy matrix Tq(u) and its transfer matrix t{u) with N 
arbitrary inhomogeneities qi 



T (u) 
t{u) = 



2u 



2u + 2 
2u+ 1 



Rw(u + <ji) • ■ • Rno(u + Qn)Rnq(u -q N )... R w {u - q x ) 



tr T (u) 



(B.5) 



2u + 2 

acting on the quantum and auxiliary spaces, and the quantum spaces, respectively. In 
addition to the YBE one has the 'reflection equation' 

i?i2 (u - v)Ti [u)R\2 (u + v)T 2 (v) — T 2 {v)Ru (u + v)T x (u)Ri 2 (u — v) 



R^u - v)T2(u)R^(u + v)Tf(v) - T»i$(u + v)T2(u)R°f(u - v). (B.6) 



The commutativity of the transfer matrices [t(u),t(u)] = can be verified by taking 
the trace of the reflection equation over the auxiliary space. When applying the 
Algebraic Bethe Ansatz it is convenient to consider the monodromy matrix as a matrix 
in auxiliary space, with entries acting as operators on the quantum spaces, namely 

/ n( u ) t}(u) 

T (u) = 



T?(«) Ti(u) 



\ Tf(u) Tf(u) 

Define the vacuum state |VAC) = |l)yi® . . .<8> |l)yiv, the state in which all sites are in 
the same state (for convenience we label the internal states by 1, . . . , N'). The action 
of the monodromy matrix T (v) on this vacuum state is 

/ a(«)|vac) i:::> ... |— ) \ 

B(u)|VAC) ... : 



(B.7) 



T(u)|VAC) = 



V 



o 







fl(u)|VAC) / 



Periodic and hard wall fermions and bosons 



24 



with the scalar functions A(u) and B(u) evaluated below. Note that entries are zero, 
except for the main diagonal (T? has the vacuum state as eigenvector, with vacuum 
eigenvalues A(u) and B(u)) and the operators Tj on the first line, which create some 
unknown state. The entries in the first line are creation operators Bj{u) := Tj(u). 
The reflection equation leads to more complicated commutation relations j between 
elements of the monodromy matrix compared to the periodic case. 

For notational simplicity we restrict the calculation from here on to the special 
case N' = 2, corresponding to particles with at most two internal spin states. In 
this case there is only one creation operator, B{u) := B2(u). Reading off the general 
commutation relations for arbitrary N' leads to a nested form j3 1351 of the Bcthc 
Ansatz. Here [B(u),B(v)] = with 

T, (v)B(u) = —, — : tB(u)T 1 (v) 

1 w w b(u~v)a(u + v) w 1 w 

" « ' 

c(u — v)b(u + v) „, ,„i . v 1 „. . ,„ . 

-h( T \ B(.v)T}(u)~ . B{v)T*(u) (B.8) 

b[u — v)a(u + v) a(u + v) 

2f « B « = — — — S(u)T|(u) 

" v ' 

B(v)Ti(u) + , tt BWT^u) 



(u-u)(u + i; + 1) w ^ w (u - w)(u + v + 1) 
2 



v)(u + V + 1) 



B{u)Tl{v). (B.9) 



The Algebraic Bethe Ansatz provides a systematic procedure for obtaining a 
complete set of cigenstatcs and eigenvectors of the transfer matrix t(u) = Tl(u) + 
T$ (u). As there is only one creation operator B(v) for the su(2) case, the ansatz for 
the Bethe state is 

M 

|$>=J]B( Mi )|VAC) (B.10) 

which is shown to be an eigenstate of t(u) for certain complex numbers {^i, . . . , /xm} 
and all arguments u. The strategy is to commute the operators T?(u) of the main 
diagonal through to the right of the creation operators B(fj,i) to act directly on the 
vacuum state. As only terms proportional to |4>) can be part of the eigenvector, 
only terms not exchanging the arguments (highlighted in the commutation relations 
above) can contribute to the eigenvector - the 'wanted terms'. All other terms, 
arising from the remaining terms in the commutation relation in any commutation 
step, are 'unwanted terms' which can be shown to vanish for certain {fii, . . . ,/j-m}- 
These conditions arc the Bethe equations. In the periodic case only one term in 
the commutation relation is 'wanted', making it possible to collect all wanted terms 
from the commutation operations in one product. Thus a simple linear transform 
T 2 2 (u) = (2u + l)T|(u) — Tl(u) is applied to obtain the commutation relations in the 
same familiar form 

„i , > a(u — v)b(u + v)„, ,„■■/ . b(2u) „. ,„■, . , 

T i (»)*(«) = 77 ( T ,{ B{u)T}(v)+ ; h( 1 . B(v)Tl(u) 

b(u — v)a{u + v) a(2u)b(v — u) 

wanted term 

| Obtained by taking special values of the external indices and evaluating the known values of the 
ij-matrix in IB. 61 . E.g., iijkl) = (2211), (2111), (2221) for the three commutation relations given. 
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a(u + v)a(2u) 



B{v)Tj{u) 
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(B.ll) 



Ti{v)B{u) 



a(v — u)a(u + v + 1) 
b(y — u)a(u + v) 



B(u)T*(v) + 



a{2v + 1) 
a(2u)b(u — v) 



B(v)T*(u) 



wanted term 
a{2v + l)6(2u) 



B(v)Tl{u). 



(B.12) 



a(u + v)a(2u) 

Combining everything together the transfer matrix eigenvalue A(u) is given by 



r(u)|*) 



2u+ 1 
2m + 2 

Al(«) 



2u + 2 Tl(u) + -±-fl(u) 



2u + 1 
1 



2u + 2 



2u + 1" 



1$) 



(B.13) 



A(«) 

with the 'wanted' contributions from the diagonal elements 

N M 



= II a ( u - 1i) a ( u + Qi) J 



ajfii - u)b(fij + u) 
b(ni - u)a{m + u) 



Kl{u) = Q{u)\ 



vacuum eigenvalue commutation relation 
1 a(u — (Mi) a(u + fii + 1) 



L J^ b(u - fii) a(u + fii) 



(B.14) 



(B.15) 



commutation relation 
The vacuum eigenvalue Q(u) of T% (u) has the required simple form if calculated as 

N 

(2u + l)T 2 2 (u) — U a(u + qi)a(u - q, 

3=1 



f| H| VAC) = 



| VAC). 



(B.16) 



The vacuum eigenvalue Q2(u) of T 2 2 (u) is obtained via additional commutation 
relations, using the YBE at v = —it, resulting in 



Q 2 («) 



c(2u) 
a(2u) 



JV AT 

I o(u + ©)a(u - qi) - JI 6(it + qt)b(u - qi) 
i=i 3=1 



N 



+ Y\_K U + Qi)K u - Qi)- 

3 = 1 

The final expression for the eigenvalue is given by 



(B.17) 



A(«) 



2u y~t u — m + 1 u + fii + 2 t-i- w + U — Qi 

11 „ _ _|- //_■ -I- 1 11 



2U + 2 1 | 11 — yUj U + |Uj + 1 1 ^ U _|_ g. _|_ 1 y _ g. -|- i 

Z— 1 2 — 1 

M 



n 



U- fli-lu + fii 



(B.18) 

±U + fli + lU - fii 

In principle it needs to be shown that the unwanted terms vanish and that the 
conditions for this are the Bethe equations. For simplicity we content ourselves 
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here with reading off the conditions which make the eigenvalue A(it) pole-free at the 
parameters /i,-, yielding the Bcthc equations in the form (with i = 1, . . . , M) 

N 

Mi 



M M 



Mi + Mj 



J- IM - fXj + 1 V: Hi + /ij + 



n 



= -i. 



(B.19) 



Applying the shifts /i& — > iAfe-c 1 — | and 
the ^-interaction literature, namely 



Mi + <Zj Mi - '/., 

ifc;c _1 gives the form normally used in 



M 

n 



Aj - A, + ic Aj 



A; - A, - ic A, 



A, 



A, 



ic 



n 



A 4 - % - |ic Ai 



A, 



|ic A, 



(B.20) 



The Bethe equations (|B.20ll are the condition for 1B.18|) to be the eigenvalue of 
the transfer matrix t(u) defined in l]B.5fl . Note that in order to make contact with 
the eigenvalue formulas in |Appendix A.2| one needs to apply the shift u — ► u — i in 
the above formula. 

The eigenvalue A HW (it) for the general su(N) case can be obtained in similar 
fashion, with 



A HW (u) 



M 

n 



(o) 



1 



(0) 



,(0) 



,(0) 



N 

£ 



2u- 1 



2u 



2u + k-22u + k-l 



M k - 

>< n 



(fe-l) | i -\ (k- 

' + k — 1 u — /X] 



1 M k 
1 -, r M 



n 



/if } + fc - 2 u - /if } - 1 



(fe) 
,(*) 



+ fc- 1 



(fe) 



(B.21) 



This formula is adapted from Ref. (SHI, with fj,f :> — > qpiA^ c 1 — (fc — l)/2, with 
the upper (lower) sign for bosons (fermions) and the monodromy inhomogeneities 
qj = /ij°' + i. Again all higher terms in the eigenvalue (i.e., sum terms with k > 1) 



vanish for the special argument u = /if and 



/' 



(0) 



1, due to terms of the form 



Y[f=i( u ~l J,[ j' > )( u + Mj U; + !)■ Here the conditions for the unwanted terms to cancel, or 
equivalently, for the eigenvalue to be pole-free, are given by 



,(oh 



(0) 



A/fc-i ^(fc) 



n 



A 



(fe-i) 



ic A 



(fe) 



A 



(fc-i) 



M k A (fc) 



A 



(fe) 



A 



(fe-i) 



M, 



a: 
(fe) 



(fc-i) 



n tw 

7=1 A l 



A 



(fe+i) 



n 



A 



CO 



A 



(fe) 



A 



(fe) 



A 



(fc) 



A 



(fe) 



±ic A 



(fe) 



A 



(fc+i) 



ic A 



(fe) 



A 



(fe) 



ic 



ic 



A 



(fc+i) 



±ic A (fc) 



A 



(fe+i) 



(B.22) 



with k = 1, . . . , N' — 1 and i = 1, . . . , Affe in the fc-th equation. Here the first and last 
Bcthc equations follow from this compact form on setting Mq = N and Mjv — 0, as 



there are no Bcthc roots A 



(JV) 



of level AJ'. 
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